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Abstract 


Symplectic unitary representations for the Poincare group are studied. The formalism is based on the noncommutative structure of 
the star-product, and using group theory approach as a guide, a consistent physical theory in phase space is constructed. The state 
of a quantum mechanics system is described by a quasi-probability amplitude that is in association with the Wigner function. As 
a result, the Klein-Gordon and Dirac equations are derived in phase space. As an application, we study the Dirac equation with 
electromagnetic interaction in phase space. 
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1. Introduction 

The first formalism to quantum mechanics in phase space 
was proposed by Wigner in 1932 [1]. He was motived by the 
problem of finding a way to improve the quantum statistical 
mechanics, based on the desity matrix, to treat the transport 
equations for superfluids [2-4]. Since then, the formalism pro¬ 
posed by Wigner has been applied in different contexts, such as 
quantum optics [5, 6], condensed matter [7-9], quantum com¬ 
puting [10-12], quantum tomography [13], plasma physics[14— 
19]. Wigner introduced his formalism by using a kind of Fourier 
transform of the density matrix, p(q, q'), giving rise to what in 
nowadays called the Wigner function, fw(q, p), where (q, p) are 
coordinates of a phase space manifold (T). The Wigner func¬ 
tion is identified as a quasi-probability density in the sense that 
fw(q, p) is real but not positive defined, and as such cannot 
be interpreted as a probability. However, the integrals cr(q) = 
f fw(q, P)dp and cr(p) = f f w (q,p)dq are distribution func¬ 
tions [1, 2]. 

In Wigner formalism each quantum operator A in the Hilbert 
space is associated with a function a w (q, p) defined in T. The 
application D.w ■ A —> aw(q,p) is such that associative algebra 
of operators in 77 defines an associative but noncommutative 
algebra in T. The noncommutativity stems from nature of the 
product between two operators in 77. Given two operators A 
and B, we have the mapping Q : AB —» aw(q, p) ★ bw(q, p), 
where the star (or Moyal)-product ★ is defined by [20] 

1 d d d d 

a w {q,p ) * b w (q, p) = a w (q,p)exp[-( — -- -z~-z~)]b w (q, p). 

2 oq op op oq 
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(Throught this Letter we use natural units: h — c — 1). Note 
that Eq.(l) can be seen as an operator A = a w (q, /;)* acting 
on functions bw(q,p), such that A{bw ) - aw * bw■ In this 
sense, we can study unitary representations of Lie groups in 
phase space using the Moyal product as defined by the operators 
A. This gives rise, for instance, to the Klein-Gordon and Dirac 
equations written in phase space [21-24]. The connection with 
Wigner function is derived, providing a physical interpretation 
for the representation. As a consequence, these symplectic rep¬ 
resentations are a a way to consider the Wigner methods on the 
bases of symmetry groups. In the present work, we apply this 
symplectic formalism to solve Dirac equation with electromag¬ 
netic interaction in phase space. These results provide a starting 
point for our analysis of nonclassical electromagnetic radiation 
sates in phase space. 

The presentation of this Letter is organized in the follow¬ 
ing way. In section 2, we define a Hilbert space 77(r) over a 
phase space T with its natural relativistic symplectic struture. 
In section 3, we study the Poincare algebra in 77(T) and the 
representation for spin 1/2. In section 4, the Dirac equation 
in phase space with electromagnetic radiation is considered. 
Quasi-amplitudes of probabilities are derived . In section 5, 
final concluding remarks are presented. 

2. Hilbert Space and Symplectic Structure 

Consider M an n-dimensional analytical manifold where 
each point is specified by Minkowski coordinates qf 1 , with p — 
0,1,2,3,4 and metric specified by diag(g) = (—h ++). The 
coordinates of each point in T*M will be denoted by (</', p >r ). 
The space T*M is equipped with a symplectic struture by intro¬ 
ducing a 2-form 

to = dq M A dp^, (2) 
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called the symplectic form (sum over repeated indices is as- and 

sumed). Consider the following bidifferential operator on C°°(T*M): W 2 = W 11 W u , (14) 



_ LJL 

(3) 


<V dPp dp 1 dq p ’ 

such that for C°° 

functions, / = /(<// p^) and g 

= g(q M ,P M ), we 

have 

{/.*} = tj(fA,gA)=fAg, 

(4) 

where 

if ] - d f d 8 d f d S 

’ 8 <V dp p dpt* dq p ' 

(5) 


is the Poisson bracket and /A and g A are two vector fields given 
by liA = Xh = —{h, ). The space T* endowed with this symplec¬ 
tic structure is called the phase space, and will be denoted by 

r. 

The notion of Hilbert space associated with the phase space 
T is introduced by considering the set of square integrable func¬ 
tions, </>in T, such that 

dpf 1 duf<p*((f,pf 1 ')^>(ci l ,pf 1 )<oo. (6) 

Then, we can write <p(cf, p M ) = (qp^\(f>), with 

J'dp^dq^iq^.p^Xq 11 ,^] = 1, (7) 

to be (0| the dual vector of \<p). We call this the Hilbert space 

WV). 


where P 2 stands for the mass shell condition and W 2 for the 
spin. 

To determine the Klein-Gordon field equation, we consider 
a scalar representation in 77(0. In this case, we can use the 
invariant P 2 to write 

P 2 < P(q^ ’P^ 1 ) = (p 2 ) * <P(qP[P), 

= (jP * p,i*)<!>{q t \p i ‘), 

= m 2 cf>(q^, fP), (15) 

where m is a constant fixing the representation and interpreted 
as mass, such that the mass shell condition is satisfied. Using 
Eq. (8), we obtain 

{ lfp » ~ V-Sr ~ ° (16) 

which is the Klein-Gordon equation in phase space. 

The association of this representation with Wigner formal¬ 
ism is given by [22] 

fw(rf l ,P >l ) = <P(q t ‘,P l ‘) * 0*(<f (17) 

where fwiq^^P 21 ) is the relativistic Wigner function. 

The representation for spin-1/2 leads to 

/ [ph - l ~ j , fP) = mi!/(<f,iP), (18) 


which is the Dirac equation in phase space, where the -/'-matrices 

3. Poincare Algebra and Dirac Equation in Phase Space fulfill the usual Clifford algebra, (y ,, y v - y v y' r ) = 2g flv . The 

„ Wigner function, in this case, is given by [22] 

Using the .star-operators, A = a w (q, p)+, we define 4-momentum 

and 4-position operators, respectively, by fP) = 'l'{(f,P >1 ) * Hq 1 ", \P ), (19) 


exp«) 


* <9) 


From Eqs. (8) and (9), we can introduce the quantity M vcr = 
Qfi P :r - QcrPfi- The operators P fI and M va are defined in the 
Hilbert space, 77(1 ), constructed with complex functions in the 
phase space E, and satisfy the set of commutation relations 


[M llv , Pa-] = i(gy,a-Pp - go-pPy), ( 10 ) 

[P ll ,Pyi= 0, (ID 

I 2d ).!V> !T j) ] — i(g iH MvcT gvp 47yyrr + gfur PI pv gva 2d p>p ) . (12) 

This is the Poincare algebra, where M pv stands for rotations and 
P fl for translations (but notice, in phase space). The Casimir 
invariants are calculated by using the Pauli-Lubanski matrices, 
W/j — \e llvpcr M^ cr P p , where e pvpcr is the Levi-Civita symbol. The 
invariants are 

(13) 


where ty((f,pA) — y°i//2(qV, pi‘), with \\P((f,pA) being the Her- 
mitian of 4/{q >l , p 1 '). We point out that the CPT theorem holds 
for non-commutative theories as showed in [25]. Therefore, 
such a theorem is also valid in phase space since the group 
structure remains the same. 

One central point to be emphasized is that the approach de¬ 
veloped here permits the calculation of Wigner functions for 
relativistic systems with methods, based on symmetry, similar 
to those used in quantum field theory. 

4. Solution of Dirac Equation with Electromagnetic Inter¬ 
action on Phase Space 

In this section, we study interactions of a spin-1/2 charged 
particle with an external electromagnetic field in Phase Space. 
The relevant equation is the Dirac equation with minimal cou¬ 
pling 

{y i P p +my F = 0, (20) 

being 

Pp -» Pfi ~ eA p , 


P 2 = P M P p , 


2 


( 21 ) 



the minimal coupling prescription, where A' = \e ljk B jX^ and 
A 0 = 0, which represents the chosen gauge. We also chose the 
magnetic field as B = (0,0, B). Thus, we have 

f/(^-eAj + m]'P = 0. (22) 

Now, we make the definition 

T = [/(/^-^)-,4// = °. (23) 


In order to obtain the energy levels, we substitute Eq. (23) into 
Eq. (22) to give 


with 

where 


\/Y ~ eA^j ( P v - eA v ) - m 2 ] /A = 0, 

yy = ^ , 


^ v = 5(/r v -rV) = ^[/,r v ]. 


The components cr 0 ', cr' 2 of the operator (26) are 


( 0 

(T ‘) 

1, (T ,] = ~ f 


0 J 

\ 


(24) 

(25) 

(26) 

(27) 


Note that these components are also expressed as cr= ia ] , 
cr'-' = -CjjijA ■ These results are explicitly evaluated in the fol¬ 
lowing representation of the y-matrices: 


1 0 
0 -1 

0 cr'' 

cr' 0 


2* = 


0 cr' 

—cr' 0 

< T k 0 

0 cr k 


75 


0 1 
1 0 


(28) 


with cr = (cr x , cr y , cr z ) being the Pauli matrices. Equation (24) 
can be written as 


[P^P M - e {P^A M + A^Pfi) + e 2 A M A II + ecP JV [P v , A /( ] - m 2 ) i/z = 0 . 

(29) 

Confining the motion at the plane XY by the choice P 3 = 0 and 
using the operators P^ = p„ - and X^ = we get 

the following equation 


d 1 d 2 


dt 4 dt 2 


~ E -iE— + - — -m \lj/ + ip x + p - — \ —r + —r 


1 ( d 2 


d 2 


eB 


if d d \ 1 I 

- I Py~ - Px~ - I + T I 


dp 


dPy 


y 4 \dx 2 dy 2 
d 2 d 2 


, d d , 


4\dydp x dxdp y 
[ xpy - yp 


\ 1 I d d 

\%- y K 


e 2 B 2 


1 d 
x H—-— 

2 dp. 


i d 

+ I y +- 

P 2 8py 


+ ieBcr 12 = 0. (30) 


Since 75 commutes through all terms of Eq. (30) and, if we 
have found the solution 1 //, then we must also have that y=,ik is a 
solution. In this case. In this case, the wave function can takes 
on one of the forms 



' <Ai ' 



# = 

^2 
^3 

. * = 

1A2 

-1A1 


, 1A4 - 


. -^2 / 


Thus, we can select only one of these solutions as that will make 
the other redundant. We specialize the solution 


y 5 iA = i/c (32) 

We can write if/ in terms of T in Eq. (23) as 

la + ys)'?, (33) 

where I is the unit matrix. From Eq. (31), we can show that Eq. 
(33) can be put in the form 


<p{E,t)4>(px,Py,X,y ) 
-ip(E,t)<p(p x ,p v ,x,y) y 
(34) 

whereyfE, t, p x , p y , x, y) is a two-component wavefunction. Note 
that, in the representation (34), the upper two components of 
Eq. (30) are now completely decoupled from the lower two. 
So, we have, in two-component form, the following equations: 


X(E,t, p x ,p y ,x,y ) 
~X(E,t,p x ,p y ,x,y) 



1 d 2 
4dP 


m 1 (p 


A 2 tp, 


(35) 


1 / d 2 


d 2 

y 4 \ dx 2 dy 2 
eB\{xp y -yp x ) + -lp y 


M Py 


d 
dp x 

+ ^ 

21 dy y dxJ 4 \dydp x 


dy 1 
d 

Pxy— 

dpy 

d 2 

dxdpy 


dx 


e 2 B 2 


x + 


2 dp. 


y + 


2 dp. 


+ ieBcr 


12 


4 2 <A, 


(36) 


where A is a constant. We point out that E is not associated to 
z'4 a priori in Eq. 35, thus the energy comes from A. Note that 
Eqs. (35) and (36) determines x and, hence, from Eq. (33), it 
determines <//. Performing a changing of variables in Eq. (36) 
of the form 


2 2 

: Px + Py + 


e 2 B 2 


eB ( yp x - xp y ) + —— (x 2 + y 2 ) , 


we note that the imaginary part of the equation vanished which 
yields 

z<t> - e 2 B 2 0 - e 2 B 2 z<p = (a 2 + seB ) 0, (37) 

where we have used icr U (p = with s = +1. If we use 
a) = zl eB and (!) = exp (-w)T’(w), the equation for F(a>) is 
found to be 


o>f" + (1 - 2oS)F - (1 - k)F = 0, (38) 

where F’ = and k = (A 2 + seB)/eB. Equation (38) is of the 
confluent hypergeometric equation type 

zF”(z) + (b-z)F'(z)-aF(z) = 0. (39) 


3 


























In this manner, the general solution for Eq. (38) is given by 

where M(a, b, z), U(a, b, z.) are the Kummer functions, and A m , 
B,„ are constants. Since only U (a,b,z) is square integrable, we 
consider it as a physical solution. Thus, we can impose that 
A m = 0. Furthermore, if a - -n, with n - 0, 1,2,..., the series 
U(a, b, z) becomes a polynomial in z of degree not exceeding n. 
From this condition, we can write 

l-k = -2n, (41) 

from which, we can extract the relation 

A 2 = eB (2 n + 1 + s). (42) 

The wave function is given by 

fm(z) = B m m(-w, 1, (43) 

where B,„ is a normalization constant. 

The Wigner function related to Dirac equation with an elec¬ 
tromagnetic interaction is formally given by 

fw(x, y, p x , py) = T„(x, y, p x , p y ) * ¥„(■*, y, p x , p y ). 

Thus Wigner function is used to determine mean values, for 
example, this result can be useful for theoretical and applied ar¬ 
eas, such as: quantum optics, quantum tomography and quan¬ 
tum computing. We point out that Fandau levels which appear 
in expression (42) represent as a matter of fact a planar oscilla¬ 
tor and the variable z in Eq. (43) give us information about the 
symplectic structure. 

5. Conclusion 

We have set forth a symplectic representation of the Poincare 
group, which yields quantum theories in phase space. We have 
derived the Klein-Gordon and Dirac equations in phase space 
and, as illustrations, studied the Dirac equation with electro¬ 
magnetic interaction. The symplectic representation is con¬ 
structed on the basis of the Moyal or star product, an ingredient 
of noncommutative geometry. A Hilbert space is then defined 
from a manifold with the features of phase space. The states are 
represented by a quasi-amplitude of probability, a wave func¬ 
tion in phase space, the definition of which makes connection 
with the Wigner function, i.e., the quasi-probability density. 
Nontrivial, yet consistent, the association with the Wigner func¬ 
tion provides a physical interpretation of the theory. Analo¬ 
gous interpretations are not found in other studies of represen¬ 
tations in phase space [25, 26]. One aspect of the procedure 
deserves emphasis. Our formalism explores unitary represen¬ 
tations to calculate Wigner functions. This constitutes an im¬ 
portant advantage over the more traditional constructions of the 
Wigner method, which entail several intricacies associated with 
the Liouville-von Neumann equation. Furthermore, the formal¬ 
ism we have described opens new perspectives for applications 
of the Wigner function method in quantum field theory. This 
aspect of the formalism will be discussed in a forthcoming pa¬ 
per. 
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